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Abstract 

The ranges of arctic-alpine species have shifted extensively with Pleistocene cli- 
mate changes and glaciations. Using sequence data from the trnH-psbA and trnT- 
trnL chloroplast DNA spacer regions, we investigated the phylogeography of the 
widespread, ancient (>3 million years) arctic-alpine plant Oxyria digyna (Polyg- 
onaceae). We identified 45 haplotypes and six highly divergent major lineages; 
estimated ages of these lineages (time to most recent common ancestor, Tmrca) 
ranged from ~0.5 to 2.5 million years. One lineage is widespread in the arctic, a 
second is restricted to the southern Rocky Mountains of the western United States, 
and a third was found only in the Himalayan and Altai regions of Asia. Three other 
lineages are widespread in western North America, where they overlap extensively. 
The high genetic diversity and the presence of divergent major cpDNA lineages 
within Oxyria digyna reflect its age and suggest that it was widespread during much 
of its history. The distributions of individual lineages indicate repeated spread of 
Oxyria digyna through North America over multiple glacial cycles. During the Last 
Glacial Maximum it persisted in multiple refugia in western North America, in- 
cluding Beringia, south of the continental ice, and within the northern limits of 
the Cordilleran ice sheet. Our data contribute to a growing body of evidence that 
arctic-alpine species have migrated from different source regions over multiple 
glacial cycles and that cryptic refugia contributed to persistence through the Last 
Glacial Maximum. 



Introduction 

Changing climates during the last 2 million years have pro- 
foundly influenced the distribution and abundance of tem- 
perate, boreal, and arctic species. Many species have under- 
gone major range shifts in response to changing sea levels or 
glacial advances and retreats. During the Last Glacial Max- 
imum (LGM) approximately 18,000 years ago, the Fenno- 
Scandian Ice Sheet in northern Europe and the Laurentide 
and Cordilleran Ice Sheets in North America extended over 
large areas of these continents (Ehlers and Gibbard 2007). 
Many tundra plants nevertheless have broad present-day dis- 
tributions in both arctic and alpine environments, a pattern 
which may reflect repeated cycles of expansion, contraction, 
and migration over several glacial cycles. Molecular studies in 
combination with fossil evidence have shown diverse phylo- 



geographic patterns among these widespread species, demon- 
strating how they came to occupy their present ranges and 
where refugia occurred (Brochmann et al. 2003; Brochmann 
and Brysting 2008). Both arctic refugia and southern moun- 
tain ranges have been suggested as source regions for arctic 
plants (Murray 1995; Brochmann and Brysting 2008). How- 
ever, the extent and timelines of contributions to the arctic 
flora from these different source areas are still not well un- 
derstood. 

Beringia (Yukon to eastern Siberia) provided a large high- 
latitude ice-free area during the Pleistocene and has long been 
proposed as a major refugium for many arctic-alpine plant 
species during glacial periods (Hulten 1937). This hypothesis 
is well supported by molecular evidence (e.g., Tremblay and 
Schoen 1999; Abbott et al. 2000; Abbott and Comes 2003; 
Alsos et al. 2005; Eidesen et al. 2007; Skrede et al. 2009). 
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Although Hulten (1937) considered Beringia the most impor- 
tant refugial region, he suggested that the broad distributions 
of many tundra species implied dispersal from additional 
refugia. Many arctic-alpine species persisted in refugia south 
of continental ice in both Europe and North America (Abbott 
and Brochmann 2003; Brochmann et al. 2003; Birks 2008; 
Brochmann and Brysting 2008); examples are Arabis alpina 
(Koch et al. 2006) and Dryas octopetala (Skrede et al. 2006). 
Additional refugia may also have existed along glacial mar- 
gins, in areas such as coastal eastern Greenland (Funder 1979) 
and coastal British Columbia (Heusser 1989; Ogilvie 1989; 
Hebda et al. 1997; King et al. 2009). Recent genetic evidence 
has indicated survival of tundra plant and animal species in 
unexpected refugia within the overall margins of continental 
ice sheets of the Late Pleistocene, both in Europe (Wester- 
gaard et al. 20 1 1 ) and in western North America (Loehr et al. 
2006; Marr et al. 2008; Shafer et al. 2010a). For tundra plant 
species, refugia have been suggested in northwestern Europe 
(Westergaard et al. 2011), northeastern Canada (Tremblay 
and Schoen 1999), northwestern Canada (Eriksen and Topel 
2006), and northern interior British Columbia (Marr et al. 
2008), but further work is needed to determine their overall 
importance. 

Genetic evidence of rapid spread has been found in some 
arctic-alpine plants. For example, the ranges of Vaccinium 
uliginosum (Alsos et al. 2005; Eidesen et al. 2007) and Rubus 
chamaemorus (Ehrich et al. 2008) expanded rapidly following 
the LGM. Long-distance dispersal and colonization is sup- 
ported by molecular evidence for a number of arctic species 
(Brochmann and Brysting 2008), and is the most plausi- 
ble explanation for current genetic characteristics of many 
plants found on the remote archipelago of Svalbard (Alsos 
et al. 2007). However, the potential for rapid migration does 
not preclude the possibility of survival in multiple locations, 
including cryptic refugia previously unknown from fossil or 
palynological evidence (Stewart and Lister 2001; Provan and 
Bennett 2008) . For example, molecular evidence suggests that 
some European tree species had a more complex glacial his- 
tory than previously thought, including survival through the 
LGM in unexpected locations (Birks and Willis 2008). Tun- 
dra plants also had complex responses to past climate change 
(Schonswetter et al. 2005), and may have spread multiple 
times from steppe or montane habitats (Tkach et al. 2008). 
An understanding of how dispersal and refugia shaped their 
present distributions can provide keys to assessing their po- 
tential responses to future climate change. 

Oxyria digyna Hill is a widespread arctic-alpine peren- 
nial herb (Fig. 1) that occurs in tundra habitats over much 
of the northern hemisphere (Fig. 2), including nearly all of 
arctic North America and Eurasia and mountain ranges of 
southern Russia, the Himalayan region, southern Europe, 
and western North America (Hulten 1971). Its fossil record 
extends back to the late Tertiary, with macrofossils dated at 




Figure 1. Oxyria digyna, Mt. Lassen, California. 



~3 million years reported from arctic Canada (Ellesmere, 
Meighen, and Prince Patrick Islands, Beaufort Formation; 
Mathews and Ovenden 1990). It grows in a variety of mesic 
tundra environments, including scree with little or no soil 
development. The flowers of Oxyria digyna are inconspic- 
uous and wind-pollinated. The winged seeds are dispersed 
moderately well by wind (Tackenberg and Stocklin 2008) and 
the species establishes effectively on glacial forelands (Whit- 
taker 1993; Stocklin and Baumler 1996). Arctic and alpine 
populations show ecotypic differences in flower and rhizome 
production and in growth responses to temperature and day 
length (Mooney and Billings 1961; Heide 2005). Oxyria dig- 
yna is taxonomically distinct, with no closely related taxa, and 
is generally treated as a single well-defined species. Oxyria 
digyna is largely diploid (2n = 14), with a few hexaploids 
(2w = 42) reported from Russia (Fedorov 1974; Goldblatt 
and Johnson 1979-2010). 

The history of the tundra floras of the northern hemisphere 
has been much more extensively investigated in Europe than 
in North America. In particular, phylogeographic studies of 
arctic-alpine plants with a focus on western North Amer- 
ica are still relatively few (Jorgensen et al. 2003; Bain and 
Golden 2005), and the contributions of refugia hypothesized 
for this region (especially cryptic refugia) are not well un- 
derstood. In a previous study (Marr et al. 2008), we used 
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Figure 2. Worldwide geographic distribution of Oxyria digyna (green 
outline) and maximum extent of ice sheets during the Last Glacial Max- 
imum (blue shading). Range map is redrawn from Hulten (1971) and 
glacial boundaries from Ehlers and Gibbard (2007). 

restriction fragment analysis to identify patterns of genetic 
diversity in western North American populations of Oxyria 
digyna in relation to recent glacial history and the locations 
of refugia. Here we used cpDNA sequence data to reconstruct 
the broader northern hemisphere phylogeographic history of 
Oxyria digyna, taking advantage of the much greater power 
of sequence data to resolve phylogenetic relationships among 
intraspecific lineages and reveal their early history. We used a 
combination of analytic approaches together with increased 
haplotype resolution and greater sampling around the north- 
ern hemisphere to address the following questions: (1) what 
is the geographic structure of genetic variation in Oxyria dig- 
yna and how does it compare with other arctic species? (2) 
How did the species and its major lineages diversify over time 
and disperse into their present geographic range? (3) What 
was the contribution of different refugia to the persistence 
of different lineages in western North America during the 
LGM? 

Methods 
Sampling 

We sequenced 181 individuals of Oxyria digyna from 140 lo- 
calities in North America, Europe, and Asia (Table 1). Sample 
localities were concentrated in western North America, but 
we obtained samples from much of the geographic range. 
We usually sampled one plant per population, but included 
1-2 additional individuals in some populations from a pre- 
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vious study that used different methods (Marr et al. 2008). 
We also sampled four outgroup taxa: Oxyria sinensis Hemsl., 
Rheum palmatum L., Rumex lapponicus Czernov, and Rumex 
paucifolius Nutt. (Table 1). 

DNA extraction, amplification, and 
sequencing 

DNA was extracted from silica-dried leaves (or in a few cases, 
pressed specimens) from 1-3 plants per population, using 
either a modified CTAB method (Doyle and Doyle 1987) 
or DNeasy Plant Mini Kits (Qiagen, Valencia, CA). We am- 
plified and sequenced two cpDNA noncoding regions, the 
trnH-psbA and trnT-trnL spacers, in all samples including 
outgroups. We also sequenced a third region, the ndhJ-trnF 
spacer, in a subset of 37 Oxyria digyna samples and in the out- 
groups. For the trnH-psbA region we used the primers psbA 
(Sang et al. 1997) and trnH (Tate and Simpson 2003). For the 
trnT-trnL region we used forward primers Tab A (Taberlet et 
al. 1991) or trnA2 (Cronn et al. 2002) with reverse primers 
Tab B (Taberlet et al. 1991) or Oxy88 (Marr et al. 2008). For 
the ndhJ-trnF region we used forward primer ndhj (Shaw et 
al. 2007) and reverse primer trnF (Dumolin-Lapegue et al. 
1997). 

PCR was performed in 50 /xL volumes with the following 
reagents: 5 /xL template DNA (1:10 dilutions), 5 /xL lOx PCR 
Buffer, 5 fiL 2 mM dNTPs, 0.25 /j,L Taq DNA polymerase 
(New England Biolabs, Ipswich, MA), and 2.5 /xL of 5 /xM of 
each primer. All PCR reactions were carried out using Techne 
TC-312 thermal cyclers. PCR parameters for amplification of 
the trnH-psbA and trnT-trnL regions were 3 min at 94°C, 30 
cycles of 94°C for 30 s, 57°C for 1 min, and 72°C for 1 min, 
with a final 10 min at 72°C. Parameters for the ndhJ-trnF 
region were 5 min at 80°C, 30 cycles of 95°C for 1 min, 50°C 
for 1 min, and 65°C for 4 min, with a final 5 min at 65°C. 
In preparation for sequencing, PCR products were purified 
using the QIAquick PCR Purification Kit (Qiagen, Valencia, 
CA). All samples were sequenced in both primer directions 
using the amplifying primers. DNA sequencing was carried 
out by Macrogen Inc. (Seoul, Korea), using ABI 3700 or ABI 
3730XL PRISM capillary sequencers (Applied Biosystems, 
Foster City, CA). To minimize PCR errors, all rare sequence 
variants were reamplified and resequenced. 

Data analyses 

Sequences for each DNA region were aligned using ClustalX 
(Thompson et al. 1997) and Jalview (Waterhouse et al. 2009), 
and manually adjusted to minimize the number of length 
variations. Three of the four outgroup taxa were highly di- 
vergent from Oxyria digyna in all sequenced regions and were 
not used further. Oxyria sinensis aligned satisfactorily with 
Oxyria digyna for the ndhJ-trnF region only. We performed 
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Table 1. Collection locality information for all samples and combined trnH-psbA and trnT-trnL haplotypes for Oxyria digyna. Localities are arranged 
by continent and geographic region, and within these in order of decreasing latitude. 

Country, locality Population ID Latitude Longitude Collector Year 1 Clade and haplotype 

Oxyria digyna 
Europe and Asia 



China, Yunnan, Gaoligong 


YN 


25 


13' N 


98° 28' E 


P. Fritsch 


2006 


F-42 


Shan 
















Greenland, Johannes V. 


GN 


83 


24' N 


29° 39' W 


D. Mazzuchi 


2006 


A-7 


Jensens Land 
















Greenland, Kap KObenhavn 


KAP 


82 


24' N 


22° 12' W 


D. Mazzuchi 


2006 


A-7 


Greenland, Zackenberg 


ZAC 


74 


28' N 


20° 34' W 


K. Westergaard 


2007 


A-7 


Station 
















Greenland, Osterdalen Valley 


OST 


69 


15' N 


53° 31 'W 


K. Westergaard 


2006 


A-14 


Greenland, Narsaq 


NAR 


60 


57' N 


46° 03' W 


K. Westergaard 


2006 


A-14 


Iceland, Skalafellsjokull 


SKA 


64 


13' N 


1 5° 44' W 


M. Reynolds 


2011 


A-14 


Nepal, Annapurna Sanctuary 


A 


28 


32' N 


83° 55' E 


W. Mackenzie 


2005 


F-41 


Nepal, Khumbu Himal, Mt. 


NEP 


27 


55' N 


86° 45' E 


A. C. Byers 


1984 


F-43 


Everest National Park 
















Norway, Endalen Valley, Sval- 


END 


78 


11' N 


1 5° 46' E 


K. Westergaard 


2006 


A-14 


bard 
















Norway, Stranda, Mt. Dal- 


NO 


62 


03' N 


7° 15' E 


O. Finch 


2003 


A-14 


snibba 
















Russia, Altai Republic, Rudnik 


RUD 


50 


20' N 


87° 45' E 


K. L. Marr, R. J. Hebda 


2010 


F-45 I 


Mine/Korumduairya Creek 
















Russia, Altai Republic, Karakjel 


KRK 


50 


02' N 


88° 02' E 


K. L. Marr, R. J. Hebda 


2010 


F-44 


Lake/Karasu Creek 
















Russia, Altai Republic, Talaya 


TAL 


50 


00' N 


89° 1 6' E 


K. L. Marr, R. J. Hebda 


2010 


F-45 


Pass 
















Russia, Kamchatka, Chazma 


KM 


58 


01' N 


1 58" 39' E 


W. Merilees 


2004 


A-4 


River 
















Russia, Magadan Oblast, near 


POB 


64 


28' N 


144" 55' E 


K. L. Marr, R. J. Hebda 


2011 


A-l 


Pobieda 
















Russia, Magadan Oblast, di- 


IKP 


63 


27' N 


146° 27' E 


K. L. Marr, R. J. Hebda 


2011 


A-l 


vide between Indygirka and 
















Kolyma Rivers 
















Russia, Magadan Oblast, 


BUR 


62 


42' N 


1 48" 48' E 


K. L. Marr, R. J. Hebda 


2011 


A-l 


Burkhalinsky Pass 
















Russia, Magadan Oblast, Ta- 


TAY 


61 


10' N 


152" 00' E 


K. L. Marr, R. J. Hebda 


2011 


A-l 


laya Pass 
















Russia, Magadan Oblast, Ar- 


ARM 


59 


40' N 


1 50° 37' E 


K. L. Marr, R. J. Hebda 


2011 


A-8 


mansky Pass 
















Russia, Sakha Republic, near 


UST 


64 


29' N 


143" 15' E 


K. L. Marr, R. J. Hebda 


2011 


A-l 


Ust Nera 
















Russia, Wrangel Island 


WRA 


71° 


15' N 


1 79" 40' W 


A. Benedict 


2009 


A-9 


Russia, Yamalo-Nenetsky Au- 


YA 


65 


27' N 


66° 50' E 


I. G. Alsos 


2004 


A-12 


tonomous 
















Region, Chernaya Mts. 
















North America 
















Canada, Nunavut, Alexandra 


AF 


78 


58' N 


75° 45' W 


H. Guest 


2003 


A-6 


Fjord, Ellesmere 
















Canada, Nunavut, Ellef 


RG 


78 


47' N 


103" 32' W 


M. Reynolds 


2005 


A-7 


Ringnes Island 
















Canada, Nunavut, Pond Inlet 


PON 


72° 


45' N 


77° 35' W 


P. Smith 


2011 


A-7 


Canada, Nunavut, Roche Bay 


RY 


68 


29' N 


82° 39' W 


B. Apland 


2006 


A-7 


Canada, Nunavut, High Lake 


HL 


67 


23' N 


110° 51' W 


B. Apland 


2006 


A-7 


Camp 
















Canada, Nunavut, Apex, Baf- 


BN 


63 


44' N 


68° 26' W 


L. Gillespie 


2004 


A-7 



fin Island 
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Country, locality 



Population ID Latitude Longitude 



Collector 



Year 



'Clade and haplotype 



Canada, NWT, Prince Patrick 
Island 

Canada, NWT, Banks Island, 

DeSalis Bay 
Canada, NWT, Palmer Lake 

Site 9 

Canada, NWT, Dechen la' 

North Canol Road 
Canada, YT, Shingle Point 
Canada, YT, Iwavik National 

Park 

Canada, YT, Upper Fish Creek, 

Richardson Mts. 
Canada, YT, Mt. Klotz camp 
Canada, YT, Ogilvie Mts. (site 

3) 

Canada, YT, Bonnet Plume, 
Gillespie Lake 

Canada, YT, Tombstone Terri- 
torial Park, North Fork Pass 

Canada, YT, MacMillan Pass 

Canada, YT, Asi Keyi Special 
Management Area 

Canada, YT, Outpost Mt. 

Canada, YT, Kotaneelee 
Range 

Canada, BC, Tagish Highland, 

TeePee Mt. 
Canada, BC, Liard Plateau, 

Caribou Mts. 
Canada, BC, Ruby Mt. 
Canada, BC, Chuck Creek 
Canada, BC, Horseranch 

Range 

Canada, BC, Cassiar Mts., 

"Dakota" Lake 
Canada, BC, Kawdy Plateau, 

High Tuya Lake 
Canada, BC, Muskwa Ranges, 

Terminal Range 
Canada, BC, Nimbus Mt. 
Canada, BC, Muskwa Ranges, 

Nonda Creek 
Canada, BC, Stikine Ranges, 

Four Mile Creek 
Canada, BC, French Range, 

Rath Mt. 
Canada, BC, Stikine Ranges, 

Little Blue Sheep Lake 
Canada, BC, Dome Mt. 
Canada, BC, Snow Peak 
Canada, BC, Three Sisters 

Range 

Canada, BC, Stikine Ranges, 

Shea Mt. 
Canada, BC, mountains south 

of Pitman River 



U L 


7C= 1 A' M 

/ D 1 H In 


1 1 0/ 

i i y 


?n' \a/ 
zu VV 


t\A Pawn i-i I He 
1 VI . rAfcfyi l(JI(Jb 




ZUU J 


A 7 


DS 


71° 27' N 


121' 


52' W 


L. O II Itripifcr 




2004 


A-7 


PLK 


64° 32' N 


129 


34' W 


A. Veitch 




2007 


A-6 


DEC 


63° 24' N 


129' 


36' W 


B. Krieckhaus 




2007 


A- 5 


SH 


69° 59' N 


137' 


25' W 


B. Bennett 




2005 


A-8 


IV 


69° 37' N 


1 40' 


50' W 


R Ronnott 
D. Del II IfcrlA 




2005 


A-8 


UF 


67° 54' N 


136' 


33' W 


B. Bennett 




2006 


A-8 


KLO 


65° 22' N 


140' 


11' W 


B. Bennett 




2007 


A-6 


G3 


64° 47' N 


137' 


43' W 


H. Guest 




2005 


A-3 (2) 


BT 


64° 44' N 


133' 


59' W 


B. Bennett 




2005 


A-l 


TB 


64° 35' N 


138' 


17' W 


K. L. Marr 




2004 


A-6, A-8, 


MP 


63° 17' N 


130' 


09' W 


B. Apland 




2006 


A-3 


AK 


61° 34' N 


140' 


50' W 


B. Bennett 




2004 


A-6 


OM 


60° 59' N 


1 38' 


25' W 


H Guest 




2004 


A-6 


KO 


60° 13' N 


124 


07' W 


B. Bennett 




2004 


A-l 


TEE 


59° 44' N 


134 


40' W 


K. L. Marr, R. J 


Hebda 


2008 


A-2 


CRM 


59° 43' N 


125° 


28' W 


K. L. Marr, R. J 


Hebda 


2007 


A-l 


RB 


59° 42' N 


133' 


22' W 


K. L. Marr 




2004 


A-6, B-15 


CC 


59° 42' N 


136' 


36' W 


H. Guest 




2004 


A-l I 


HR 


59° 28' N 


1 28' 


56' W 


IS.. L. 1 VI d 1 1, r\. J 


Hebda 


2004 


A-6 


UL 


J3 ZD IM 


1 3rv 

I JU 


1 3' \A/ 
I J VV 


IS.. L. IVIdll, r\. J 


ntuud 


ZUU4 


A 1 A in D 1C 
H I, M- 1 U, D 1 D 


HT 


59° 1 5' N 


130' 


3 1 ' W 


k" 1 KAxrr R 1 
IS. L. Ivldl 1, r\. J 


Hebda 


2004 


A-l B-15 


TER 


59° 02' N 


125' 


55' W 


IS. L. Ivldl 1, r\. J 


Hebda 


2007 


A-l 


NIM 


59° 02' N 


132' 


30' W 


K. L. Marr, R. J 


Hebda 


2008 


B-15 


N 


59° 00' N 


125' 


30' W 


K. L. Marr 




2003 


A-1, B-18 (2) 


4MC 


58° 51' N 


129 


3V W 


K. L. Marr, R. J 


Hebda 


2007 


A-6 


RTM 


58°45' N 


130' 


25' W 


K. L. Marr, R. J 


Hebda 


2007 


A-6 


LBS 


58° 44' N 


128' 


15' W 


K. L. Marr, R. J 


Hebda 


2007 


A-6 


DM 


58°31' N 


129' 


35' W 


K. L. Marr, R. J 


Hebda 


2004 


A-6 (2) 


SP 


58° 20' N 


130' 


23' W 


W. Mackenzie 




2003 


B-15 


SI 


58° 20' N 


129' 


45' W 


W. Mackenzie 




2003 


A-6 


SHM 


58° 1 9' N 


128' 


57' W 


K. L. Marr, R. J 


Hebda 


2007 


A-6 


PIT 


57° 53' N 


127' 


54' W 


K. L. Marr, R. J 


Hebda 


2009 


A-6 
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Country, locality 


Population ID 


Latitude 


Longitude 


Collector 




Year 


'Clade and haplotype 


Canada, BC, Thudaka Range 


THU 


57° 48' N 


126° 


37' 


W 


K. L. Marr, R. J. 


Hebda 


2009 


B-15 


Canada, BC, Spatsizi Provincial 


SZ 


57° 41 ' N 


128° 


52' 


W 


R. J. Hebda 




2005 


A-6, B-15 


Park, Airplane Creek 




















Canada, BC, Mt. Edziza 


E 


57° 38' N 


130° 


40' 


W 


C. Howers 




2005 


B-15 (2), B-17 


Canada, BC, Moosehorn Lake 


MOO 


57° 35' N 


127° 


13' 


W 


K. L. Marr, R. J. 


Hebda 


2009 


B-15 


Canada, BC, Brothers Lake 


BRO 


57° 12' N 


127° 


26' 


W 


K. L. Marr, R. J. 


Hebda 


2009 


E-34 


Canada, BC, Kemess North 


K 


57° 04' N 


126° 


45' 


W 


J. Pojar, W. Mackenzie 


2003 


A-l 


Canada, BC, Wrede Range 


WRE 


56° 38' N 


126° 


12' 


w 


K. L. Marr, R. J. 


Hebda 


2009 


B-15 


Canada, BC, Oweegee Range 


ow 


56° 34' N 


129° 


28' 


w 


K. L. Marr, R. J. 


Hebda 


2006 


B-17 


Canada, BC, Chase Mt. 


C 


56° 33' N 


125° 


16' 


W 


K. L. Marr, R. J. 


Hebda 


2003 


B-15 


Canada, BC, Mt. Henri 


MT 


56° 29' N 


124° 


45' 


W 


K. L. Marr, R. J. 


Hebda 


2003 


E-34 


Canada, BC, Ludington Peak 


LP 


56° 28' N 


123° 


23' 


w 


K. L. Marr, R. J. 


Hebda 


2003 


B-18 


Canada, BC, Needham Creek 


NE 


56° 24' N 


123° 


31' 


w 


K. L. Marr, R. J. 


Hebda 


2003 


A-l, B-15, B-18 


Canada, BC, Hanna Ridge 


HN 


56° 14' N 


129" 


26' 


w 


K. L. Marr, R. J. 


Hebda 


2006 


B-17 


Canada, BC, Mt. Tommy Jack 


TO 


56°03' N 


127° 


46' 


W 


K. L. Marr, R. J. 


Hebda 


2006 


E-34 


C snsrl; 3 RC Mnrnhpp Mt 

^ d I OUd, UV., IVIUI J C~ VIL. 


MR 


55=26' N 


123° 


02' 


W 


K. L. Marr, R. J. 


Hebda 


2003 


E-40 


Canada, BC, Murray Mt. 


M 


55° 24' N 


122° 


37' 


w 


K. L. Marr, R. J. 


Hebda 


2003 


E-39 


Canada, BC, Mt. Crum, "Holz- 


H 


55° 01' N 


121° 


31' 


w 


K. L. Marr, R. J. 


Hebda 


2003 


E-40 


worth Meadows" 




















Canada, BC, Babine Mts. 


BA 


55° 00' N 


126° 


56' 


W 


K. L. Marr, R. J. 


Hebda 


2002 


E-34 (2) 


Canada, BC, Nass Ranges, Mt. 


CU 


54° 54' N 


128° 


42' 


w 


K. L. Marr, R. J. 


Hebda 


2006 


E-34 


Couture 




















Canada, BC, Sleeping Beauty 
Mt. 


SB 


54° 34' N 


128° 


52' 


w 


M. Cheney 




2003 


B-17 (2) 


Canada, BC, Mt. Thornhill 


TL 


54° 31' N 


1 28" 


27' 


W 


G. Allen 




2004 


B-17, E-34 


Canada, BC, Nadina Mt. 


NA 


54° 06' N 


126° 


52' 


W 


K. L. Marr, R. J. 


Hebda 


2002 


E-34 (2) 


Canada, BC, Tweedsmuir Peak 


TW 


53° 39' N 


126" 


29' 


w 


K. L. Marr, R. J. 


Hebda 


2002 


D-32 (2) 


Canada, BC, Canoe Mt. 


CA 


52° 43' N 


1 19" 


1 1 ' 


w 


K. L. Marr 




2002 


E-35 (2) 


Canada, BC, Mt. De La 


LT 


52° 42' N 


132" 


02' 


W 


G. W. Douglas 




2003 


B-17 


Touche. 




















Canada, BC, Itcha llgachuz 


CH 


52° 42' N 


1 24° 


51 ' 


W 


R. Coffey 




2003 


E-35 


Canada, BC, Mackenzie Mt. 


MK 


52° 35' N 


126° 


09' 


w 


W. Mackenzie 




2003 


E-35 


Canada, BC, Tweedsmuir Park, 


CL 


52° 32' N 


125° 


49' 


w 


R. Coffey 




2003 


D-32, E-35 


Crystal Lake trail 




















Canada, BC, Caput Mt. 


CAP 


52° 25' N 


120° 


36' 


W 


F. Osorio 




2006 


E-40 


Canada, BC, Wells-Gray 


TR 


51° 47' N 


119° 


55' 


W 


K. L. Marr 




2002 


E-40 


Provincial Park, Trophy Mt. 




















Canada, BC, Niut Range 


NT 


51° 43' N 


124° 


39' 


w 


W. Mackenzie 




2003 


E-35 (2) 


Canada, BC, Brisco Range, 


KIN 


50° 42' N 


116° 


00' 


w 


K. L. Marr 




2007 


E-34 


Kindersley Pass 




















Canada, BC, Bugaboo Pass 


BO 


50°41 ' N 


1 16° 


45' 


w 


O. McDadi 




2004 


E-34, E-40 


Canada, BC, Paradise Mine 


PA 


50° 29' N 


1 16° 


18' 


w 


K. L. Marr 




2002 


D-32, E-34 


Canada, BC, Vancouver Island, 


MW 


50° 21' N 


127" 


18' 


w 


K. L. Marr, R. J. 


Hebda 


2005 


E-34 


Merry Widow Mt. 




















Canada, BC, Kootenay Lake 


KOO 


50° 1 8' N 


117° 


08' 


W 


W. Mackenzie 




2007 


E-39 


Canada, BC, Whistler Mt. 


WH 


50° 04' N 


122° 


57' 


w 


K. L. Marr 




2003 


E-35 


Canada, BC, Manning Park 


MA 


49° 6' N 


120° 


46' 


w 


K. L. Marr 




2002 


E-40 (2) 


Canada, BC, Gimli Peak 


Gl 


49° 46' N 


117" 


39' 


w 


H. Guest 




2006 


E-35 


Canada, BC, Vancouver Island, 


FP 


49° 41 ' N 


125" 


22' 


W 


K. L. Marr 




2003 


E-34 


dUUVfcr l_l[Llfcrl LdNtr 




















Canada, BC, Fernie, Three Sis- 


3S 


49° 35' N 


115" 


08' 


w 


K. L. Marr 




2002 


E-35 (2) 


ters trail 




















Canada, BC, Vancouver Island, 


AR 


49° 14' N 


124" 


36' 


W 


K. L. Marr 




2004 


D-32 


Mt. Arrowsmith 




















Canada, BC, Silver Tip Mt. 


SM 


49° 1 0' N 


121" 


13' 


W 


G. W. Douglas 




2003 


E-35 


Canada, BC, Kishinena Range 


KR 


49° 02' N 


114" 


16' 


w 


W. Mackenzie 




2005 


E-35 


Canada, AB, Pipestone R. Trail 


PP 


51° 35' N 


116" 


09' 


w 


S. White 




2003 


E-35 
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Table 1. Continued. 



Country, locality 



Population ID 



Latitude 



Longitude 



Collector 



Year 



'Clade and haplotype 



LF 


51" 31 


N 


ii6° irw 


S. White 


2003 


E-35 


BFS 


69° 58 


N 


144° 46' W 


J. Jorgensen 


2007 


A-13 


ANA 


68" 08 


N 


151° 36' W 


B. Krieckhaus 


2008 


A-6 


DK 


63"39' 


N 


144° 37' W 


H. Guest 


2005 


A-10 


RH 


63-21' 


N 


145° 42' W 


H. Guest 


2005 


A-l (2) 


CR 


62° 37 


N 


143° 28' W 


H. Guest 


2005 


A-10 


CM 


61 "34' 


N 


144" 36' W 


H. Guest 


2005 


A-3 (2) 


CW 


61" 02 


N 


149- 07' W 


H. Guest 


2005 


A-l 


HSM 


60" 33 


N 


172° 42' W 


CAN 587881 


1944 


A-l 


CO 


60° 27 


N 


143° 53' W 


A Batten 


2005 


B-18, B-19 


KF 


59° 27 


N 


149° 57' W 


M. Carlson 


2005 


B-15 


CF 


57° 42 


N 


135° 02' W 


B. Krieckhaus 


2006 


B-15 (2) 


KD 


57° 16 


N 


154° 12' W 


C. Parker 


2005 


B-15 


UNF 


53° 51 


N 


166° 32' W 


B. Krieckhaus 


2007 


B-16 


WIN 


48° 57 


N 


121° 39' W 


K. L. Marr 


2007 


E-34 


BAG 


48" 51 


N 


121° 41' W 


K. L. Marr 


2007 


E-34 


HP 


48" 44 


N 


1 20° 40' W 


K. L. Marr 


2003 


E-35, E-40 


OP 


47° 55 


N 


123° 23' W 


K. L. Marr 


2003 


B-21, D-32 


1 R 


46" 52 


N 


1 2 1 " 3 1 ' W 


G. Allen, K. L. Marr 


2004 


E-39 


SY 


46" 03 


N 


1 13° 58' W 


H. Guest 


2006 


E-35 


AB 


45" 40 


N 


1 1 0° 34' W 


P. Lesica 


2004 


B-23, B-25 


HY 


45" 26 


N 


110" 57' W 


K. L. Marr 


2005 


D-30, E-37 


HG 


45" 22 


N 


1 1 6° 30' W 


G. Allen, K. L. Marr 


2004 


B-24 


SD 


45°21' 


N 


1 1 6" 3 1 ' W 


G. Allen, K. L. Marr 


2004 


E-36 


EG 


45" 24 


N 


121° 40' W 


G. Allen, K. L. Marr 


2004 


B-20 (2) 


TJ 


44° 24 


N 


121° 48' W 


G. Allen, J. Antos 


2004 


E-39 


SW 


44° 18 


N 


11 8° 41' W 


G. Allen, K. L. Marr 


2004 


B-22, E-38 


AS 


42" 05 


N 


122° 43' W 


G. Allen, J. Antos 


2004 


D-31 


HM 


42" 44 


N 


109" 10' W 


J. Dunn 


2006 


C-28 


MD 


41" 21 


N 


106" 20' W 


H. Guest 


2006 


C-27 


FL 


39" 55 


N 


105° 40' W 


K. L. Marr 


2003 


C-26, C-29 


L 


40" 29 


N 


121° 31' W 


G. Allen, J. Antos 


2004 


D-31 


PU 


37" 14 


N 


118° 39' N 


H. Guest 


2006 


E-35 




57" 59 


N 


98° 27' E 


L. Q. Jun 




2006 










G. Allen 




2008 




56° 49 


N 


128- 35' E 


K. L. Marr, R. J. Hebda 




2005 




43° 17 


N 


1 09° 50' E 


G. Allen, J. Antos 




2005 



Canada, AB, Lower Fish Lake 
USA, AK, Beaufort Sea coast, 

Camden Bay 
USA, AK, Anaktuvuk Pass, 

Alaska 
USA, AK, Dot Lake 
USA, AK, Richardson Highway 
USA, AK, Caribou Creek 
USA, AK, Chitina Rd. 
USA, AK, Crow Pass 
USA, AK, Bering Sea, St. 

Matthew Island 
USA, AK, Carbon Mt. 
USA, AK, Kenai Fjords 
USA, AK, Chichagof Island, 

Hoonah Ridge 
USA, AK, Kodiak Island 
USA, AK, Unalaska Island, 

Pyramid Peak 
USA, WA, Winchester Mt. 

Lookout 
USA, WA, Mt. Baker, Bagley 

Lakes Trail 
USA, WA, Hart's Pass 
USA, WA, Olympic Mts., Ob- 
struction Point 
USA, WA, Chinook Pass 
USA, MT, Anaconda Pintlar 

Range, Seymour Lake 
USA, MT, Absaroka Mts. 
USA, MT, Hyalite Creek 
USA, ID, Heaven's Gate 
USA, ID, Seven Devils Mts. 
USA, OR, Mt. Hood, Elliott 

Glacier 

USA, OR, Three-Fingered Jack 
USA, OR, Strawberry Mts. 
USA, OR, Mt. Ashland 
USA, WY, Wind River Range, 

Haystack Mt. 
USA, WY, Medicine Bow Mts. 
USA, CO, Forest Lakes 
USA, CA, Mt. Lassen 
USA, CA, sw of Bishop, Paiute 

Pass Trail 
Outgroups 
Oxyria sinensis 
China, Yunnan, Deqing 
Rheum palmatum 
Canada, Jardin botanique de 

Montreal (cultivated) 
Rumex lapponicus 
Canada, BC, Anthony Creek 
Rumex paucifolius 
USA, Wyoming, Wind River 

Mts., Green River Lake 



'Parenthese indicate multiple samples from the same locality. 
© 2012 The Authors. Published by Blackwell Publishing Ltd. 
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maximum parsimony (MP) and maximum likelihood (ML) 
analyses of the ndhJ-trnF region with Oxyria sinensis as 
outgroup; the resulting trees showed little resolution of rela- 
tionships among Oxyria digyna haplotypes and are not re- 
ported here. All reported analyses are based on the trnH-psbA 
and trnT-trnL spacer regions. 

Haplotype networks for the trnH-psbA and trnT-trnL non- 
coding regions, separately and in combination, were con- 
structed using the statistical parsimony software TCS v. 1.21 
(Clement et al. 2000). An ILD test (Michewich and Farris 
1981), as implemented in PAUP 4.0M0, was carried out us- 
ing 500 random replicates to check for partition heterogeneity 
of the two regions. Indels were treated as unit characters and 
included in all analyses. To obtain estimates of support for 
branches of the combined network, we performed MP and 
Bayesian phylogenetic analyses on a data set consisting of one 
sequence for each combined haplotype (with indels coded as 
0/1 characters in a separate data partition). MP analysis was 
performed with PAUP 4.0 version blO (Swofford 1998), with 
clade support estimated from 500 bootstrap replicates us- 
ing a heuristic search with TBR branch swapping. Bayesian 
analysis was carried out with MrBayes v. 3.1.2 (Huelsenbeck 
and Ronquist 2001) under the assumption of a strict molec- 
ular clock, using a GTR substitution model as determined 
with jModelTest 0.1.1 (Guindon and Gascuel 2003; Posada 
2008). We executed two sets of four Markov chain Monte 
Carlo (MCMC) runs for 2 million generations, sampling ev- 
ery 100th generation, with burnin of 25%. 

We estimated time to most recent common ancestor 
(Tmrca) for major lineages (with >5 samples) identified in 
the phylogenetic analyses, using BEAST software (Drum- 
mond and Rambaut 2007). Analyses were based on the com- 
plete data set of 181 sequences and all characters except indels 
(which were excluded because they are difficult to incorpo- 
rate into coalescent analyses) . Two independent analyses were 
each run for 1 0 s generations under a GTR substitution model, 
selected using jModeltest 0.1.1. Analyses were combined and 
evaluated using Tracer 1.5 (Rambaut and Drummond 2007). 
In preliminary analyses assuming a relaxed (uncorrelated log- 
normal) molecular clock or exponential population growth, 
marginal posterior distributions of (respectively) clock rate 
variation and population growth rate included zero; we there- 
fore assumed a strict molecular clock and constant popula- 
tion size for all subsequent runs. In order to estimate T M rca 
in years, replicate analyses were carried out with calibration 
based on either fossil evidence or a published mutation rate. 
For fossil calibration we assigned an age of 3.0-3.4 million 
years, the estimated age of the oldest known Oxyria digyna 
macrofossil (Matthews and Ovenden 1990), to the root of the 
tree (using a normally distributed prior of 3.2 ± 0. 14 million 
years, to give a 95% confidence interval of approximately the 
above range). For calibration from mutation rate, we used 
a mean substitution rate of ~1.52 ± 0.06 x 10 -9 substitu- 
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tions per site per year as estimated for chloroplast noncoding 
regions by Yamane et al. (2006). 

We investigated possible expansion events in the histories 
of major lineages of Oxyria digyna by performing mismatch 
analyses using Arlequin v. 3.5 (Excoffier et al. 2005) on all 
lineages with five or more samples (Clades A, A + C, B, D, 
E, and F). Mismatch analyses test the hypothesis that a mis- 
match distribution (the distribution of all possible pairwise 
differences for a group of sequences) approximates the uni- 
modal shape predicted for a population or lineage that has 
undergone a demographic or spatial expansion event (Ex- 
coffier 2004); the mean number of pairwise differences can 
be used to infer time since an expansion event. 

In western North America, the geographic area for which 
we had the highest sampling density, we assessed spatial pat- 
terns of genetic diversity for five different geographic regions. 
These were delineated to approximately reflect differences in 
tundra environment (high arctic to southern alpine) and 
glacial history (glaciated or not), as follows: (1) Alaska and 
the Yukon; (2) Northwest Territory and Nunavut; (3) north- 
ern British Columbia (latitudes 54.5-60.0°N); (4) southern 
British Columbia (latitudes 49.0°-54.5°N); and (5) conti- 
nental western USA (south of 49.0°N). We tabulated overall 
numbers of haplotypes and numbers of haplotypes unique to 
a region (across all clades and within major clades) and calcu- 
lated genetic diversity, H (the probability that two randomly 
chosen haplotypes in a sample are different) and molecular 
diversity, n (the mean number of pairwise differences be- 
tween haplotypes in a sample), both within clades and overall 
clades for each region. Both diversity measures are sensitive 
to numbers and relative abundances of haplotypes; n also 
reflects divergence among haplotypes. All genetic diversity 
measures were calculated using Arlequin v. 3.5 (Excoffier 
etal. 2005). 

Results 

Sequence variation 

Total aligned length of the trnH-psbA region was 315 bp, 
with 24 variable characters (18 substitutions, five indels of 5 
to 8 bp in length, and one single -base T repeat). Total aligned 
length of the trnT-trnL region was 910 bp, with 68 variable 
characters used for analyses (41 substitutions and 27 indels 
of 3 to 19 bp); five additional 1- or 2-bp repeats could not 
be aligned unambiguously and were omitted from all anal- 
yses. There was no evidence for heterogeneity of the two 
regions (ILD test in PAUP*, P = 0.11) and all inferences are 
based on the combined two-region data set with an aligned 
length of 1225 bp. Sequences representing all distinct vari- 
ants found in this study are deposited in GenBank; accession 
numbers are JQ080918-JQ080962 for the trnH-psbA region 
and JQ080963-JQ081007 for the trnT-trnL region. 
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Phylogenetic relationships among 
haplotypes 

TCS analysis of the combined trnH-psbA and trnT-trnL data 
set yielded a total of 45 haplotypes (Fig. 3) , which formed dis- 
tinct groups corresponding with strongly supported clades 
(all with Bayesian posterior probabilities of 0.99-1.00) iden- 
tified in the phylogenetic analyses (Fig. 3a). These included 
a widespread clade (A) comprising nearly all arctic samples; 
four clades (B to E) with variously overlapping geographic 
distributions in western North America; and a clade ( F) found 
only in the mountains of central and southern Asia. A sin- 
gle, additional haplotype (26) was of indeterminate phyloge- 
netic position but had greatest sequence similarity to Clade C 
(Fig. 3b). Haplotypes in sampled populations are indicated 
in Table 1. 

Although each of the six major clades was strongly sup- 
ported, relationships among them were mostly unresolved 
in the Bayesian phylogenetic analysis (Fig. 3a). Clades A and 
C formed a well-supported higher-level clade, A + C (with 
posterior probability of 0.96; Fig. 3a). 

Timelines and clade histories 

Estimated T M rcaS for the analyzed major clades suggest that 
at least some of these have a very long history spanning 
much of the Pleistocene (Table 2). Estimated coalescence 
times were similar whether calibrated with fossil evidence 
or mutation rate. Clade F was identified as the oldest clade 
(mean T M rca 2.58-3.02 x 10 6 years). The combined Clades 
A + C (mean T M rca 1.42-1.66 x 10 6 years) and Clade D 
(meanT M RCA 1.27-1.45 x 10 6 year) appear to be substantially 
older than the remaining clades (A, B, and E), of which the 
youngest, Clade E, had a mean T M rca of 458-516 x 10 3 years 
(Table 2). Data for Clade C were insufficient for analysis. 

Mismatch analyses of major clades (A, A + C, B, D, E, and 
F) suggested spatial expansion events in the histories of all 
clades (P > 0.05 for all tests of fit to the model). Mismatch 
distributions were approximately unimodal for all of these 
clades, but with a second, small peak (suggesting the pos- 
sibility of multiple expansion events) in Clade A + C and 
in Clade D. Estimates of r had very broad 95% confidence 
limits for all clades and were insufficiently accurate for useful 
estimates of times since expansion; these are not reported 
here. 

Haplotype diversity and geographic 
distribution 

Haplotype diversity of Oxyria digyna in western North Amer- 
ica (Table 3) varied among geographic regions, both overall 
and within clades. Overall haplotype diversity measures H 
and 7r were high in Regions 1 (Alaska/Yukon), 3 (northern 
British Columbia), and 5 (continental western USA) and 
lowest in Region 2 (Nunavut and Northwest Territories). 
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Haplotype diversity was low in Region 4 (southern British 
Columbia), in marked contrast to Region 3 (northern British 
Columbia). The high overall values of it in several regions, 
especially Regions 3 and 5, reflect the presence of haplotypes 
belonging to divergent major clades in this region (Table 3). 
Haplotypes from several clades were present in all regions ex- 
cept Region 2. Within individual clades, haplotype diversity 
varied with latitude across geographic regions, tending to be 
highest either in the north (Region 1) or in the south (Region 
5) (Table 3). These two geographic regions also harbored 
unique (region-specific) haplotypes for all major clades that 
were present, and unique haplotypes were concentrated in 
these regions (Table 3). 

Three of the major lineages appear to have very disparate 
geographic distributions, with little overlap (Fig. 4a) . The sin- 
gle Eurasian lineage, Clade F, is genetically and geographically 
remote from all other haplotype groups detected (although 
it was undoubtedly undersampled in this study). Clade C 
was found only in a small region of the Rocky Mountains 
in the western USA. The grouping of this small clade with 
the northern Clade A in the phylogenetic analysis suggests 
a formerly more contiguous distribution over the interven- 
ing region. The largely arctic Clade A is the most widely 
distributed lineage (Fig. 4a, b), spanning a continuous longi- 
tudinal range of at least 195° (from eastern Greenland across 
North America to eastern Siberia) , and extending to 270° with 
inclusion of the Yamalo-Nenetsky region of Russia. Clade A 
comprises a subclade of closely related haplotypes (1 to 13) 
and a single divergent haplotype (14) of trans-Atlantic dis- 
tribution (Greenland, Svalbard, and northern Europe). The 
main subclade is widely distributed, with a concentration of 
haplotypes in northwestern North America. Within this sub- 
clade, the earliest diverging haplotypes 1 and 13 were found 
from southwestern Alaska to northern British Columbia 
(Fig. 3b, 4b), with the highest diversity of haplotypes also 
in this region. 

Clades B, D, and E all occur in western North America, 
mostly south of Alaska and the Yukon, with broad geograph- 
ical overlap (Table 3, Fig. 4a, b). Their relative abundances 
and the distribution patterns of haplotypes within each clade 
suggest different histories. Clade B includes 11 haplotypes 
divided into two well-marked and geographically distinct 
subclades (Fig. 3, Fig. 4b). The predominantly northern sub- 
clade (haplotypes 15 to 21) occurs primarily from south- 
ern coastal Alaska to northern British Columbia, with a dis- 
junct southern lineage (haplotypes 20 and 21). A second, 
southern subclade (haplotypes 22 to 25) is restricted to cool, 
damp alpine habitats in the western USA. Clade D is a sparse 
clade with four divergent haplotypes (30 to 33) distributed 
over a large geographic region in western North America 
(Fig. 4b). Haplotype 33 in this clade, identified from a single 
northern Yukon locality, is sister to a clade of the remaining 
three haplotypes, found further south. Clade E with seven 
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Yukon, Alaska, E Russia 
W Russia 

NE Canada, Greenland 
Greenland, Europe 
N BC 
S BC 
W USA 
Central Asia 



Figure 3. Relationships among Oxyria digyna 
haplotypes from the combined analysis of 
trnH - psbA and trnT-trnL cpDNA spacer 
regions. (A) Bayesian phylogenetic tree, with 
clade support indicated by Bayesian posterior 
probabilities and maximum parsimony 
bootstrap values. (B) TCS haplotype network. 
Circle diameters are proportional to the 
number of samples of each haplotype; 
connecting dots indicate inferred haplotypes 
not found in the data set. Colors indicate the 
geographic occurrences of haplotypes. 
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Table 2. Estimates of time to most recent common ancestor (T MRCA ) in years before present for major lineages of Oxyria digyna, based on combined 
trn H- psb A and trn T-L cpDNA sequence data. 





Tmrca (from fossil date) 1 


Tmrca (from substitution rate) 2 


Lineage 


Mean ± SD 


Mean ± SD 


0. digyna (all extant clades combined) 3 




3,705,000 ±821,000 


Clade A 


769,000 ± 279,000 


893,000 ± 266,000 


Clades A + C 


1,41 5,000 ± 545,000 


1,656,000 ± 583,000 


Clade B (all) 


777,000 ± 294,000 


91 1,000 ±288,000 


Clade B (northern subclade) 


475,000 ± 191,000 


561,000 ± 190,000 


Clade D 


1,265,000 ±518,000 


1,454,000 ± 507,000 


Clade E 


458,000 ± 204,000 


516,000 ± 199,000 


Clade F 


2,576,000 ± 696,000 


3,024,000 ± 1,188,000 



1 based on late Tertiary macrofossils of O. digyna (Matthews and Ovenden 1 990). 
2 based on substitution rate of Yamane et al. (2006) for cpDNA noncoding regions. 
Calculated from substitution rate only. 



haplotypes (34-40) is geographically well-defined, occurring 
from north-central British Columbia to central California. 
Clades B, D, and E overlap in a large contact zone extending 
over much of the formerly glaciated region of western North 
America. 

Discussion 

An old and complex history 

Oxyria digyna is an old and morphologically distinct species 
(Murray 1995; Brochmann and Brysting 2008). Its diversity 



of chloroplast haplotypes is consistent with its long fossil 
history extending back to the late Tertiary (Fredskild and 
R0en 1982; Mathews and Ovenden 1990; West 2000; Akopian 
et al. 2008; Bennike et al. 2010). Chloroplast haplotype di- 
versity in Oxyria digyna is greater than in most other widely 
distributed arctic-alpine species investigated to date, includ- 
ing V. uliginosum (Alsos et al. 2005; Eidesen et al. 2007), A. 
alpina (Koch et al. 2006), and Juncus biglumis (Schonswetter 
et al. 2007), and is comparable to genetic diversity of multi- 
species lineages in genera such as Hordeum (Jakob and Blat- 
tner 2006). Alsos et al. (2005) identified three major cpDNA 



Table 3. Oxyria digyna genetic diversity statistics by geographic region for localities in western North America, based on combined trn H- psb A and 
trn T-L cpDNA sequence data. Diversity measures for each region are given for all haplotypes combined (in bold) and for haplotypes within each major 
clade. 



Region 


Clade 


No. of samples 


No. of haplotypes 


No. of unique haplotypes 


Gene diversity, H 


Molecular diversity, n 


Region 1 (Alaska/Yukon) 


All 


29 


11 


6 


0.894 


4.259 




A 


23 


6 


3 


0.834 


1.652 




B 


5 


4 


2 


0.900 


2.400 




D 


1 


1 


1 


0 


0 


Region 2 (Nunavut/NWT) 


A (All) 


9 


3 


2 


0.556 


1.667 


Region 3 (northern BC) 


All 


51 


11 


2 


0.842 


6.753 




A 


21 


5 


2 


0.638 


0.895 




B 


22 


3 


0 


0.507 


1.221 




E 


8 


3 


0 


0.607 


1.643 


Region 4 (southern BC) 


All 


39 


6 


0 


0.775 


4.289 




B 


4 


1 


0 


0 


0 




D 


5 


1 


0 


0 


0 




E 


30 


4 


0 


0.653 


1.076 


Region 5 (w USA s of 49°N) 


All 


26 


20 


15 


0.979 


6.200 




B 


7 


6 


6 


0.952 


4.762 




C 


3 


3 


3 


0.909 


3.000 




D 


4 


3 


2 


0.833 


3.333 




E 


11 


7 


3 


0.909 


1.927 




hap. 26 


1 


1 


1 


0 


0 
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Figure 4. Geographic distributions of all haplotypes 
(Fig. 3, Table 1) from the combined analysis of trnH - 
psbA and trnT-trnL cpDNA spacer regions of 181 
samples of Oxyria digyna. (A) Major clades. (B) 
Individual haplotypes within major clades (haplotype 
26 is included on the Clade C map). Haplotypes 
occurring at multiple localities are indicated by circles 
and those found at only one locality by triangles. 
Colors are used only to differentiate haplotypes and 
do not carry geographic information. Haplotype 
numbers correspond to those in Fig. 3. 



lineages within the arctic shrub V. uliginosum, which also has 
a Tertiary fossil history, and concluded that the earliest diver- 
gence of these lineages probably occurred at the beginning of 
the Quaternary glaciations (> 700,000 years). Our estimated 
overall T M rca for Oxyria digyna based on published molec- 
ular substitution rates suggests an origin as long ago as the 
late Tertiary (~3.7 million years; Table 2), which is consistent 
with the fossil evidence. Murray (1995) suggested that Oxyria 
digyna is one of the oldest arctic-alpine species. 



Oxyria digyna shows a striking phylogenetic split into at 
least six deeply divergent cpDNA lineages. Although rela- 
tionships among these clades were not resolved in the phy- 
logenetic analyses, estimated times to most recent common 
ancestor (T M rca) for the different clades indicate that they 
are of substantially different ages. All of the major clades ap- 
pear to be old (~0.5-2.5 x 10 6 years), and some lineages 
(clade A + C, clade D, clade F) may have originated in the 
early Pleistocene. Molecular evidence increasingly indicates 
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(Brochmann and Brysting 2008) that many species in the 
relatively young (2-3 million year old) arctic tundra ecosys- 
tem harbor greater genetic variation than expected. How- 
ever, among arctic species for which data are available, only 
Saxifraga oppositifolia appears to have as much divergence 
among major cpDNA lineages (Abbott et al. 2000) as does 
Oxyria digyna. The presence of divergent lineages within a 
species indicates not only a long history, but sufficiently large 
populations that genetic diversity is maintained over time. Al- 
though shifts in climate through the Pleistocene may have led 
to genetic bottlenecks for many arctic-alpine species, highly 
cold-adapted species may have maintained very large ranges 
that simply shifted in response to climate cycles (Brochmann 
and Brysting 2008). This is a likely scenario for the very 
widespread and cold-tolerant Oxyria digyna. 

Pan-arctic taxa are of diverse origins. Some of the oldest 
species may have arisen in southern alpine habitats and others 
in long-persistent high-latitude regions such as Beringia (Ab- 
bott and Brochmann 2003; Brochmann and Brysting 2008). 
Among widely distributed species, multiple sources seem 
likely. Evidence concerning the geographic origins of Oxyria 
digyna is equivocal. Its nearest relatives occur in Asia; its 
only congener, Oxyria sinensis, is restricted to western China 
and the sister taxon of Oxyria is the wholly Eurasian genus 
Rheum (Frye and Kron 2003). However, the oldest known 
fossils of Oxyria digyna are in arctic Canada and in Green- 
land. Oxyria digyna may have originated in central or eastern 
Asia but diversified in North America early in its evolution. 
The high haplotype diversity of major clades in northwestern 
North America suggests a long history of Oxyria digyna in 
this region. The complex patterns of lineage distribution and 
overlap in western North America likely reflect a series of 
regional extinction and recolonization cycles associated with 
multiple climatic fluctuations and glacial episodes in the later 
Pleistocene. 

Histories of major lineages 

The differing geographic distributions of Oxyria digyna lin- 
eages imply different Pleistocene histories. The western North 
American Clade C and the Eurasian Clade F were found only 
in southern alpine areas. Oxyria digyna is disjunct in north- 
ern and southern Russia (Hulten 1971; Fig. 2) and Clade 
F may represent a southern lineage, although it is probably 
more widespread than our sampling indicates. Hulten (1937) 
and others have suggested that many widely distributed arc- 
tic species migrated alternately south and north as glaciers 
advanced and retreated, sometimes resulting in disjunct ge- 
ographic ranges with isolated southern populations. 

In contrast, the widespread arctic distribution of Clade A 
suggests a high latitude origin. Clades A and C form a higher- 
level clade, which could have spread from Asia or possibly 
originated in western North America with subsequent diver- 
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gence of Clades A and C to the north and south. The high 
haplotype diversity of Clade A in northwestern North Amer- 
ica and the distribution of its ancestral haplotypes are consis- 
tent with a long presence in the Beringian region, with more 
recent eastward and westward spread through arctic regions. 
Oxyria digyna is moderately well dispersed by wind (Tack- 
enberg and Stocklin 2008) and probably was able to spread 
rapidly into new habitat as continental ice sheets receded 
and climatic conditions became more favorable. Central and 
eastern regions of arctic Canada and northern Greenland are 
occupied almost exclusively by a single derived haplotype, 
suggesting recent migration. The history of arctic lineages of 
Oxyria digyna in northern Asia is less clear. Fossils indicate 
that it occurred in northern Siberia during and preceding 
the LGM (Kienast et al. 2001; Lopatina and Zanina 2006), 
and clade A haplotypes in northwestern Siberia are closely 
related to those from western North America. Broad geo- 
graphic regions of low haplotype diversity, with higher di- 
versity in refugial areas such as Beringia, have been observed 
in other arctic plants, including A. alpina (Koch et al. 2006), 
Ranunculus glacialis (Schonswetter et al. 2003), and Draba 
species (Skrede et al. 2009). The single amphi-Atlantic hap- 
lotype 14 in Clade A may represent a distinct lineage that 
persisted through the LGM in either European or northeast 
North American refugia. 

Three lineages of Oxyria digyna (Clades B, D, and E) have 
an unusually broad zone of overlap extending over much 
of British Columbia and adjacent regions. The sympatry of 
these divergent cpDNA lineages and the geographic distribu- 
tions of haplotypes within each lineage indicate a long and 
complex history of Oxyria digyna in western North America, 
involving repeated expansion, migration, and local extinc- 
tion. Mismatch analyses suggested spatial expansion events 
in the histories of all these clades, though we could not ascer- 
tain whether these events coincided with initial clade diver- 
gence or occurred subsequently. The differences in estimated 
age and relative abundance of these clades suggest succes- 
sive episodes of migration from one or more source regions. 
Clade D is an old (~1.2— 1.5 x 10 6 years), sparse clade and 
probably represents an early-colonizing lineage of which only 
remnant haplotypes now remain. Clades B and E, in contrast, 
are of more recent origin (~0.4-0.9 x 10 6 years). Clades B 
and D both show evidence of range fragmentation resulting 
from glacial advances. Clade B has a particularly complex pat- 
tern of north-south disjunctions between related haplotypes, 
suggesting multiple fragmentation events in its history. 

Persistence through the last glacial 
maximum in western North America 

In species with a long history, early distributions of lineages 
are difficult to infer because of the strong influence of later 
events. Present-day spatial patterns of genetic variation are 
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often best understood in relation to the most recent glacia- 
tion and associated climatic changes. This is true for Oxyria 
digyna, although the time since the LGM (18,000 years) rep- 
resents less than 1% of the history of the species (>3 mil- 
lion years) and thus almost all of the observed haplotypes 
must predate the LGM. Most tundra and boreal species now 
present in western North America persisted through the last 
glacial maximum north or south of the continental ice sheets 
(Jorgensen et al. 2003; Beatty and Provan 2010); some also 
persisted in refugia between these limits (Loehr et al. 2006; 
King et al. 2009). Oxyria digyna gives evidence of persistence 
in all of these areas, with different refugia harboring different 
combinations of clades. 

The ice-free area of Beringia formed a northern refugium 
extending over most of Alaska and the northern half of the 
Yukon (Fig. 2) in which Clades A and D almost certainly 
persisted through the LGM. Clade A is widespread in the 
Beringian region, which is also where its haplotype diversity 
is highest. Clade D includes a haplotype (33) that was found 
only in Beringia; the other haplotypes occur much farther 
south and likely persisted elsewhere. Clade B may also have 
persisted in Beringia or nearby cryptic refugia; it is distributed 
along the southern margin of Alaska including the Aleutian 
Islands. 

Areas south of the continental ice sheets also provided ex- 
tensive suitable habitat. Two separate lineages within Clade 
B and all haplotypes of Clade C are presently restricted to 
cool moist alpine habitats south of the limits of continen- 
tal ice. These are probably remnants of much more extensive 
populations that were widespread during the LGM and subse- 
quently became isolated into small, disjunct patches of alpine 
habitat. Clade E and the southern subclade of group D also 
occur south of the glacial boundary, but extend north into 
formerly glaciated regions of British Columbia. These likely 
persisted, at least in part, south of the continental ice sheets, 
migrating northward as the ice receded. Clade E occurs from 
the western USA north over much of British Columbia, with 
some haplotypes found almost entirely within the limits of 
the Cordilleran ice sheet but with highest diversity of haplo- 
types in the south. Northward movement from south of the 
Cordilleran ice sheet and persistence within its boundaries 
are both plausible for this clade. 

Our evidence suggests persistence of Oxyria digyna within 
the boundaries of the Cordilleran ice sheet during the LGM. 
We observed marked differences in haplotype diversity be- 
tween northern and southern British Columbia (Table 3, Fig. 
4b), although these regions, both within the boundaries of the 
Cordilleran ice sheet, are thought to have had a similar glacial 
history over the last 25,000-30,000 years. Genetic diversities 
of Clades A and B were nearly as high in northern British 
Columbia as in adjacent unglaciated regions further north, 
and included haplotypes found only within the boundaries of 
the Cordilleran ice sheet. In contrast, within-clade genetic di- 



G. A. Allen et al. 

versities in the southern part of the ice sheet are much lower, 
implying continuous ice cover at the LGM. One hypothesis 
to explain these differences is that the extensive tundra habi- 
tat of northern refugia facilitated more extensive postglacial 
migration of Oxyria digyna from the north. However, hap- 
lotype diversity of Oxyria digyna is even higher south of the 
Cordilleran ice sheet than north of it ( Table 3 ) , yet most of this 
variation is not found in southern British Columbia. Thus 
persistence of Oxyria digyna in ice-free habitats within the 
northern Cordilleran ice sheet, rather than differential sub- 
sequent migration, appears to be the most likely explanation 
for the observed pattern of haplotype diversity. Persistence in 
northern British Columbia during the LGM is also indicated 
for two vertebrates (Loehr et al. 2006, Shafer et al. 2010a), and 
may have been possible for other species (Shafer et al. 2010b). 
The occurrence of refugial habitat within the margins of con- 
tinental ice sheets in Europe is similarly supported by recent 
molecular studies (Brochmann et al. 2003; Westergaard et al. 
2011). 

Conclusions 

Oxyria digyna is a very old species, as indicated by both fossil 
evidence and the cpDNA patterns reported here. It has higher 
genetic diversity than most arctic-alpine species studied thus 
far, reflecting its age, but also suggesting a widespread distri- 
bution over time with few bottlenecks. Its history rivals that 
of multispecies complexes, reflecting both ancient and more 
recent cycles of range contraction and expansion. We show 
that different lineages within a single species may have very 
different histories and dispersal patterns spanning hundreds 
of thousands of years. Although the geographical origin of 
Oxyria digyna is most likely in Eurasia, our evidence sug- 
gests a long and complex history of this species in western 
North America. Differing ages of lineages and distributions of 
haplotypes indicate repeated migrations over multiple glacial 
cycles. Our genetic evidence indicates refugia for Oxyria dig- 
yna not only north and south of continental ice sheets but 
also within the limits of the northern Cordilleran ice sheet 
in northern British Columbia. Thus our findings support the 
emerging generalization that cryptic refugia were important 
for plant survival and recolonization of glaciated regions, as 
already demonstrated for Europe (Brochmann et al. 2003; 
Brochmann and Brysting 2008). 
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